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' We study both polarized and unpolarized proton structure functions in the kine- 

' matical region of large Bjorken a; and four-momentum transfer of few GeV^. In this 

region the phenomenon of parton-hadron duality takes place between the smooth 
continuation of the deep inelastic scattering curve and the average of the nucleon 
resonances. We present results on a perturbative-QCD analysis using all recent 
^ l' accurate data with the aim of extracting the infrared behavior of the nucleon 

^ , structure functions. 

> : 

' 1. Introduction 

^ I Parton-hadron duality is generally defined as the similarity between 

hadronic cross sections in the Deep Inelastic Scattering (DIS) region and 
in the resonance region. It encompasses therefore a range of phenomena 
where one expects to observe a transmogrification from partonic to hadronic 
degrees of freedom, a question, the latter, at the very heart of Quantum 
ChromoDynamics (QCD). 

A number of experiments were conducted in the early days of QCD 
where the onset of parton-hadron duality was observed as the equivalence 
between the continuation of the smooth curve describing different observ- 
ables from a wide variety of high energy and large momentum transfer re- 
actions - structure functions, sum rules, R{s) for e^e~ hadrons, heavy 
meson decays... - and the same observables in the low energy region char- 
acterized by low final state invariant mass values and resonance structure. 
A fully satisfactory theoretical description of this phenomenon, that be- 
came to be accepted as a "natural" feature of hadronic interactions, is still 
nowadays very difficult to obtain. Recent progress both on the theoretical 
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and experimental side has however both renovated and reinforced the 
hadronic physics community's interest in this subject as also demonstrated 
by the very lively discussions among participants in the present Workshop. 

In our contribution to the Workshop, we present evidence that stan- 
dard Perturbative QCD (PQCD) approaches in the large Bjorken x region 
might not be adequate to describe parton-hadron duality. In particular, 
by conducting an analysis of the most recent polarized and unpolarized 
inclusive electron scattering data, we unravel a discrepancy in the behav- 
ior of the extracted power corrections from the DIS and resonance regions, 
respectively. 

2. Overview of data and QCD-based Interpretation 

Parton-hadron duality being the idea that the outcome of any hard scat- 
tering process is determined by the initial scattering process among ele- 
mentary constituents - the quarks and gluons - independently from the 
hadronic phase of the reaction, is a well rooted concept in our current view 
of all high energy phenomena. Bloom and Oilman (BG) duality ^ is the 
extension of this idea to a kinematical region characterized by lower center 
of mass energies of the hard scattering process. Recently, more accurate 
experimental data have been collected that allow us to explore in detail this 
phenomenon. Besides the already mentioned new data on both inclusive 
electron-proton scattering ^, and polarized inclusive electron-proton scat- 
tering ^ , several additional data sets and hadronic reactions were measured 
in the resonance region; polarized inclusive electron-proton scattering at 
Jefferson Lab kinematics r — > v+ hadrons ^, — > 7r+n and, fi- 
nally, inclusive electron-nucleus scattering ^ (most of the recent data were 
presented and discussed at this Workshop, see also for a recent review 
of both theoretical and experimental results). 

A particularly interesting result was found in studies of inclusive reac- 
tions with no hadrons in the initial state, such as e+e^ hadrons, and 
hadronic r decays ®. It was pointed out that, because of the truncation of 
the PQCD asymptotic series, terms including quark and gluon condensates 
play an increasing role as the center of mass energy of the process decreases. 
Oscillations in the physical observables were then found to appear if the 
condensates are calculated in an instanton background. Such oscillating 
structure, calculated in ^ for values of the center of mass energy above the 
resonance region, is damped at high energy, hence warranting the onset of 
parton-hadron duality. 



February 2, 2008 4:39 Proceedings Trim Size: 9in x 6in 



fantoni 



3 

In what follows wc examine a related question, namely whether it is 
possible to extend the picture of duality explored in the higher region 
^, to the resonance region, or to the BG domain. A necessary condition 
is to determine whether the curve from the perturbative regime smoothly 
interpolates through the resonances, or whether, instead, violations of this 
correspondence occur. The latter would indicate that we are entering a 
semi-hard phase of QCD, where preconfinement effects might arise. Our 
analysis of "duality violations" requires a sufEcicutly large and accurate set 
of data. We have applied it therefore to both the unpolarized and polar- 
ized inclusive measurements of proton structure functions in the resonance 
region. 

3. Monitoring the Transition between pQCD and npQCD 

We outline two important procedures for the study of parton-hadron dual- 
ity in structure functions: i) the continuation of DIS curve into the reso- 
nance region; ii) the averaging of the resonances. Although these concepts 
are equivalently found in a number of different reactions, and in different 
channels (see e.g. ^"'^), in this contribution we concentrate on the pro- 
ton structure functions gi and F2, for polarized, and unpolarized electron 
scattering, respectively. 

3.1. Continuation of DIS Curve 

It is important to define exactly what one means by "continuation" of the 
DIS curve, in order to be able to define whether parton-duality can be con- 
sidered to be fulfilled. The accuracy of current data allows us, in fact, to 
address the question of what extrapolation from the large Q^, or asymptotic 
regime the cross sections in the resonance region should be compared to. In 
principle any extrapolation from high to low is expected to be fraught 
with theoretical uncertainties ranging from the propagation of the uncer- 
tainty on as{M^) into the resonance region to the appearance of different 
types of both perturbative and power corrections in the low regime. All 
of these aspects need therefore to be evaluated carefully. 

Our approach applies to the large Bjorken x behavior of inclusive data. 
We therefore consider: 

(a) Non-Singlet (NS) Parton Distribution Functions (PDFs) evolved at 
Next to Leading Order (NLO); 

(b) PQCD evolution using the scale « — z) which properly takes into 
account integration over the parton's transverse momentum ^^'^"^-^ 
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(c) Target Mass Corrections (TMCs). 

We perform an extensive study of inclusive data in the resonance region by 
extrapolating to this region all available parameterizations of PDFs, which 
can be considered pure DIS, down to the measured ranges for x, Q^, and 
final state invariant mass, = Q^{l—x)/x + M^. Our general approach 
for both the unpolarizcd structure function F2ix,Q'^), and the polarized 
one, gi{x,Q'^) is described in detail in An important point illustrated 
also in is that the uncertainty due to the use of different parameteriza- 
tions can be taken into account by a band that is currently smaller than 
the experimental one in the region of interest. A potential theoretical er- 
ror in the extrapolation of the ratios to low could be generated by the 
propagation of the error in as{M^). However, because at large x pertur- 
bative evolution involves only Non-Singlet (NS) distributions, we expect it 
to affect minimally the extrapolation of the initial pQCD distribution, even 
to the low values of W'^ considered. 

The problem of resumming the large logarithm terms arising at large x 
was first noticed in a pioneering paper There it was shown how this type 
of resummation can be taken care of by considering the correct definition 
of the upper limit of integration for the transverse momentum in the ladder 
diagrams defining the leading log approximation. This implies replacing 
with « W'^ — Q^{1 — z), i.e. an invariant mass, in the evolution equations. 
Such a procedure was used to obtain our results both in Refs.^^'^^, and in 
the current contribution. 

Finally, TMCs are expected to be important at W"^ M^. Although 
TMCs are kinematic in nature and their contribution to the Operator Prod- 
uct Expansion (OPE) was calculated early on their effect on structure 
functions evolution cannot be evaluated straightforwardly at leading twist 
since a truncation of the twist expansion brings inevitably to a mismatch 
between the Bjorken x supports of the TM-corrected and the "asymptotic" 
results The extent of this mismatch is, however, small in the DIS region 
thus rendering approximated treatments applicable In the resonance 
region the discrepancies arising in the ^'^ approach can be large. We adopt, 
therefore, the prescription of Ref. according to which the twist expan- 
sion for the standard moments Mellin of the structure function including 
TMCs is truncated consistently at the same order in for both the 

kinematic terms and the dynamical higher twists. As noticed already in 
this method applies for suflaciently small values of the expansion parameter 
pa 4M^a;^/Q^. In addition, we have control on the uncertainty which is a 
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term of 0(1/Q*). 



3.2. Averaging Procedure 

Resonant data can be averaged over, according to different procedures We 
considered the following complementary methods; 



where F^'^** is evaluated using the experimental data in the resonance region 
^. In Eq.(l), for each ^alue: a^min = Q^/iQ^ + W^^ax - M^), and 
aJmax = Q'^/iQ'^ + H^min - M'^). W^i^ and W^max deUmit either the whole 
resonance region, i.e. VFmin ~ 1-1 GeV^, and VF^ax ~ 4 GeV^, or smaller 
intervals within it. In Eq.(2), ^ is the Nachtmann variable and Mn{Q^) 
are Nachtmann moments The r.h.s. of Eq.(3), F^^'^^{^, W'^), is a smooth 
fit to the resonant data ^, valid for 1 < W"^ < 4 GeV^; F^'^'* symbolizes the 
average taken at the = {x, W"^) of the data. 

4. Results 

After describing our program to address quantitatively all sources of the- 
oretical errors started in i^'^^, in Fig.l we present our main results on the 
extraction of the dynamical Higher Twist (HT) terms from the resonance 
region. A clear discrepancy marking perhaps a breakdown of the twist ex- 
pansion at low values of W'^ is seen for the unpolarized structure function, 
F2 (upper panel). A comparison with other results obtained in the DIS re- 
gion ^® is also shown. For the polarized structure function, gi, we added to 
our previous analysis data from ^ at = 0.65, 1, 1.2 GeV^. In addition, we 
used the experimental values of the ratio R = (Jl/(^t from recent Jefferson 
Lab measurements in the resonance region which introduce an oscillation 
around the original result of about 2%, well within the error bars. A com- 
plete presentation and discussion of these results along with comparisons 
with other extractions will be given in a forthcoming paper Prom 



^Similar formulae hold for the polarized structure function, pi. 




(1) 





F|™(,x, Q\x, W^)) = F^'-^C, W^) 



(3) 
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the figure one can sec that although the trend seen in of a large violation 
of duality at « 1 GeV^ seems to be confirmed, more polarized data at 
large x are needed in order to draw definite conclusions. 

This work is partially supported (S.L.) by the U.S. Department of En- 
ergy grant no. DE-FG02-01ER41200. 
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Figure 1. Upper panel: Comparison of HT contributions for the structure function F2 
in the DIS and resonance regions, respectively. The full circles are the values obtained in 
the resonance region . For F2 these are compared with extractions using DIS data, from 
Lower panel: ratio of the experimental data on gi and the PQCD extrapolation 
to the resonance region ^■^ . Notice that the new results obtained using data from * agree 
with the trend of the Hermes data. 
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